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The Casimir mutual free energy F for a system of two dielectric concentric nonmag- 
04 ' netic spherical bodies is calculated, at arbitrary temperatures. Whereas F has recently 

^ , been evaluated for the special case of metals (refractive index n = oo), here analogous 

t^~- 1 results are presented for dielectrics, and shown graphically when n = 2.0. Our calcula- 

' tional method relies upon quantum statistical mechanics. The Debye expansions for the 

C D ' Riccati-Bessel functions when carried out to a high order are found to be very useful 

in practice (thereby overflow/underflow problems are easily avoided), and also to give 
accurate results even for the lowest values of /. Another virtue of the Debye expansions 
is that the limiting case of metals becomes quite amenable to an analytical treatment 
in spherical geometry. We first discuss the zero-frequency TE mode problem from a 
mathematical viewpoint and then, as physical input, invoke the actual dispersion rela- 
tions. The result of our analysis, based upon adoption of the Drude dispersion relation 
as the most correct one at low frequencies, is that the zero- frequency TE mode does not 
contribute for a metal. Accordingly, F turns out in this case to be only one half of the 
conventional value. 



^ , Keywords: Casimir effect; quantum statistical mechanics. 



1. Introduction 

Consider the free energy F{T) at temperature T due to the mutual interaction 



^ . between two spherical dielectric bodies with concentric surfaces at r = a and r = 

b. Between the two media there is a vacuum. An analysis of this problem was 
recently given in Ref. 1, using both quantum statistical methods and field theoretical 
methods. Whereas the general formalism worked out in Ref. 1 was valid for arbitrary 
values of the (equal) permittivities e in the two dielectric regions, r < a and r > b, 
the explicit evaluations of F{T) for various temperatures T and widths d = (b — a) 
in Ref. 1 were assuming perfectly conducting walls only, corresponding to e — > oo. 
Our purpose in the present paper is to extend this calculation so as to include 
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general values of the permittivity. As to our knowledge, such a calculation has not 
been undertaken before, ^Ithough there are similarities with the theory given by 
Kleinert some years ago □. We will assume, as in Rcf. 1, that the two media are 
nonmagnetic. The present account is a short version of a forthcoming paper U. 

The general conclusion drawn in Ref. 1 was that it is the quantum statistical 
method that is the most simple and powerful approach when one is to handle the 
case of general e. The most central formula in our context is the statistically 
derived Eq. (40) in Ref. 1; it gives the value of /3F = F/T for arbitrary values 
of temperature, width, and e. Whereas this equation in terms of a very compact 
notation, it will be convenient here to rewrite it slightly. Let m g (—00,00) be an 
integer corresponding to Matsubara frequencies K — 27rm//3; let n — ^ be the 
refractive index of the two media lying at r < a and r > b, and let si{x), ei{x) be 
Riccati-Bessel functions with imaginary argument defined according to 

si{x) = y^/^(x), ei{x) = \I^K^{x), (1) 

so that their Wronskian becomes W{si,ei} = —1. Here v = 1 + 1/2, and I^, K^, 
are modified Bessel functions. We write the formula as 

oc ' 00 

PF^Y. E(2^ + l)Kl-AP0+ln(l-Ar^)], (2) 

m=0 1=1 

where the prime on the summation sign means that the m = term is taken 
with half weight. The two eigenvalues Xf^ and Xf^ in Eq. (2) correspond to the 
transverse magnetic and the transverse electric modes. (In the notation of Ref. 1, 
Xei = Xf^', Xi = A^'^ .) For later use we will write these eigenvalues as ratios. 
First, 

where 

/i = ns'i{x)si(nx) - si{x)s'i{nx), 
/a = ne'i{y)ei{ny) - ei{y)e'i{ny) , 
fs = ne'i{x)si{nx) - ei{x)s'i{nx), 

fi = nei{ny)s'i{y) - e'i{ny)si{y), (4) 
a; and y being the nondimensional frequencies 

x — 2T:ma/[3, y = 2-Kmb/j3. (5) 

We put Ti = c = kg — 1. It should be emphasized that, in contradistinction to 
the formalism in Ref. 1, the primes in Eqs. (4) mean derivatives with respect to the 
whole argument. Our present way of writing is in accordance with current usage. 
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Next, the TE eigenvalues are written as 



9i92 

9m' 



(6) 



where 



91 



s'i{x)si{nx) — nsi(x)s'i{nx) 



92 = e'i{y)ei{ny) - nei{y)e'i{ny), 
53 = e'i{x)si{nx) - nei{x)s'i{nx), 
9i = ei{ny)s'i{y) - ne'i{ny)si{y). 



(7) 



The foUowing point should be noted. The formulas (4) and (7) contain, as a special 
case, the situation when the walls aX r = a,b are perfectly conducting. This case, 
as already mentioned, corresponds to setting n = ^/e — oo. When considering 
the contribution from zero Matsubara frequency, m = 0, we are confronted with a 
delicate two-limit problem. The conventional way to proceed when handling this 
problem within the framework of nondispersive theory, has been to take the limits 
in the following order: (i) Set first e = oo; (ii) take then the limit m — > 0. This 
is made analytically by observing the small-argument expressions for the Riccati- 
Bessel functions. This way of taking the limits was advocated already in the 1978 
paper of Schwinger, DeRaad, and Milton and the same procedure was followed in 
Ref. 1. If we follow the same procedure also now we get, by insertion into Eq. (2), 
the following free energy expression for perfectly conducting walls: 



m=0 1=1 



sijx) ei{y) 



e'lix) s'liy) 



(8) 



which is in agreement with Eq. (68) in Ref. 1. If we next, following the same recipe, 
let a; —> 0, ?/ — > 0, we obtain the following contribution from m = 0: 



ppconv^^ _^ OO, m = 0) = V(2Z + 1) In 



1=1 



b 



(9) 



again in agreement with Ref. 1, Eq. (79). This is the conventional result. Both the 
two electromagnetic modes are in this way found to contribute equally to the sum 
in Eq. (9). 

A discussion has recently arisen as to whether this recipe for dealing with the 
TO = term in the TE mode is really correct. The problem becomes most acute in 
the high T regime, but is present at moderate and low temperatures also. We refer 
here to the paper of Bostrom and Sernelius El, questioning this point, and the subse- 
quent comment of Lamoreaux i. What has been most welcome in recent years are 
the accurate experiments on the Casimir force, due to Lamoreaux and Mohideen 
et al. El El E3. By means of these experiments it becomes much easier to formulate 
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a sound theory. Several thei 
periments' variousJacets M 



leQie jcaj 



ers have lately appeared, discussing the ex- 



An extensive recent review has been given 
by Bordag et al. t3. We also mention several other related papers E3 Ej, of a 
more general nature; also these being concerned with finite temperature effects in 
a Casimir context. 

Below, we will make use of the present formalism to analyze how the m = case 
works out for the spherical geometry. It turns out that the formalism becomes quite 
amenable. An important ingredient in our analysis is the use of the Debye expan- 
sions for the Riccati-Bcssel functions. They make the formalism quite transparent, 
and they imply that the overflow or underflow problems that so easily turn up in 
this sort of calculations, are easily abandoned. Moreover, when carrying out the 
Debye expansion to a high order (18th order in the quantity 9 defined by Eq. (15) 
below), the numerical accuracy becomes high for all values of I, quite sufficient for 
all practical calculations. We analyze the problem of metals first from a mathemat- 
ical point of view and then, inserting the plasma dielectric model versus the Drude 
model as physical input, show how the result for F depends critically on which 
dispersion relation one chooses. As physically the Drude model is preferable at very 
low frequencies, we conclude that the m = TE mode does not contribute to F 
for a metal. This implies that the conventional expression for F for a metal has 
to be multiplied by one half. This is also in agreement with our earlier statistical 
mechanical considerations for the static case in Sec. Ill in Ref. 1. 

2. Numerical Evaluation 

2.1. The Debye expansions 

We define the nondimensional temperature: 

27ra 

p 

implying a; = mt, and write the Debye expansions of the Riccati-Bcsscl functions 
in the form EJ 

-'(-)4 [i+^]V. -^^^^^^[^(-)]' (11) 



s[{x) = 1 [1±£!M^ e-'(-) C[9{x)], (13) 
2 Vz{x) 

e[{x) = Jl±^^^ D[0{x)]. (14) 
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Here = /+ 1/2, Z = 1, 2, ... and 

z{x) = x/u, e{x) = [1 + ^2(a;)]-l/2^ 

{9 is the same as the symbol t in Ref. 20). The four polynomials, A{0), 3(9), C{9),D{e), 
are found to be of order unity. In Ref. 21 we expanded them to order 0^^. These 
expansions, which will not be reproduced here, are found to be easily handled by a 
computer. The polynomials possess the following important property: 

{A{e),B{e),C{e),D{9)} ^1 when 61^0. (16) 

The factors in Eqs. (11)-(14) that can take extreme values, are the exponentials. 
They are easily dealt with analytically in the Debye formalism. We avoid in this 
way the overflow/underflow problems that might easily occur from a simple use of 
the computer Bessel library. 

It is now convenient to calculate the following ratios between the functions 
defined in Eq. (4): 



and 



/i ^ 1 2.„rx) n^lC[9{x)] - A[9{x)]C[9{nx)]/A[einx)] 
h 2 rfi^D[9{x)]+B[9{x)]C[e{nx)]/A[e{nx)] ^ ' 

h = _2,-2..to) nHD[9{y)] - B[9{y)]D[9{ny)]l B[9{ny)] 

U nHC[9{y)]+A[6{y)]D[6{ny)]/B[9{ny)y ^ ' 



where 7 and 5 are the coefficients 



S^ll+lM. (19) 



Consequently (cf. Eq. (3)) 
where 



(20) 



V{x) - viy) = v/TT^ - VTT^) l^'^.l^'yi ] , (21) 

and where [...] is the polynomial ratio following from Eqs. (17) and (18). 
Similarly 

il = _1 -2..(.) ^C[9{^)] - A[9{x)]C[9{nx)]/A[9{nx)] 

gs 2 ^D[e{x)]+B[9{x)]C[9{nx)]/A[e{nx)]' ^ ' 

92 ^ _^-2.r,(y) ^mv)] " B[9{y)]D[9{ny)]/ B[9{ny)] 

94 5C[9{y)]+A[9{y)]D[9{ny)]/B[9{ny)]' ^ ^ 
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SO that A^'^ takes the same form as A^^^, Eq. (20), only with the difference that 
[...] now is formed from the polynomials in Eqs. (22) and (23). 

2.2. Numerical results for dielectrics 



n = 2.0 



0.5 1 1.5 2 2.5 



Fig. 1. Logarithm of nondimcnsional free energy, logjQ(— /3Ft), versus relative width d/a = [b — 
a)/a for various temperatures t = 2™//?. Refractive index n = 2.0. 

Figure 1 shows, as an example, how \ogiQ{—f3Ft) varies with logj^Q i when n = 
2.0. It shows clearly the presence of a low-temperature plateau. For higher values 
of there is a gradual change into the region where F varies linearly with t. Our 
calculations show that the magnitude \F\ of the free energy for an ordinary dielectric 
is much less than for a metal (n = oo)El. This is as we would expect. 

Generally, we found the asymptotic Debye expansions to be useful for a; > 10 
and/or I > 9. Then, an accuracy of 8 digits for the individual terms was achieved. 
Below these limits for x and I, we employed the machine-generated Bessel functions. 
For small values of d/a and t, slow convergence was observed. The summation of the 
series thus became rather demanding. As an example, when d/a = 0.05, t = 0.01, 
about 1.1 million terms were needed, if we truncated the summation at e = 10~^ 
(here e means the ratio between a general term in the series and the sum). The 
sum itself is however accurate only up to 4 or 5 digits. 

An important result was that even for low values of the asymptotic series gave 
very good results. One reason for this is the high-order expansions used for the 
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polynomials A, B, C, D. It seems most likely that the Debye expansions (at least 
when carried out to order 9^^) can be used for all x and for all practical purposes. 

3. The Limiting Case of a Metal 



3.1. The nondispersive case 

It turns OTit that the Dcbyc expansions are quite useful also for the purpose 
of analyzing metals. We first assume a nondispersive medium, so that the case of 
metals corresponds to e ^ oo, or n — > oo. 

For simplicity let us categorize how to take the two limits, i.e. the limit on n, 
and the limit on the Matsubara number m. The first option, which we shall call 
option A, means taking first n ^ oo, thereafter m — > 0. Option B reverses the 
succession of n and m. 

We consider the TM mode, first employing option A. Taking the first limit 
n oo, it follows that 0{x) is finite, whereas 9{nx) 0. Thus all polynomials 
{A, B, C, D}[6{x)] taken at argument 6{x) are finite, whereas {A, B, C, D}[6{nx)] 
1 according to Eq. (16). Since n'^-y and n^S are proportional to n for large n 
according to Eq. (19), we get from Eqs (17)-(20) 



A/^(n^oo) =e' 



Hv{^)-r,iy)] C[9{x)]D[e{y)] 
D[9{x)]C[9{y)] ■ 



(24) 



Taking the limit m ^ we have from Eq. (21) r]{x) — r]{y) ln(a/6), so that the 
m = free energy becomes 



<TM 



(m = 0) = -^(2Z + l)ln 1 - 



21+1 



(25) 



Consider next option B. When m — !■ we have 9{x) 1, 9{nx) 0, implying 
that y 1, 6 1 and \J^{m = 0) = (a/5)^'+^. Thus we obtain the same 
expression for the m = TM free energy as before, Eq. (25). The robustness of the 
TM calculated free energy is actually what we could expect on physical grounds: 
the TM mode means that the magnetic field is transverse to the radius vector r, 
thus parallel to the spherical surfaces at r = a, h. This is precisely the natural 
electromagnetic boundary condition for the magnetic field at perfect conducting 
surfaces. 

Consider then the TE mode. Employing option A we get 



AP(n ^ oo) = e2''[''(^)-''(^)l lim 



{yC[9{x)] - A[9{x)]}{5D[9{v)] B[9{v)]} 



{yD[9{x)]+ B[9{x)]}{6C[9{y)]+ A[9{y)]y 

(26) 

The difference between this case and the preceding case lies in the sensitivity of 
Eq. (26) with respect to 7 and 5. From Eq. (19) it follows that 7^0, 5^0 
implying that, when we take the limit m ^ 0, Xj^ (a/6)^'+^. Then, Eq. (2) 
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shows that the TE contribution to the to = free energy becomes the same as 
Eq. (25). 

Consider now option B. We obtain 7-^1, S ^ 1. Then, according to Eqs. (22) 
and (23), Xf^ when m 0. Consequently 

B : pF^^im 0) = 0. (27) 

Option B thus gives only half as large total free energy as the conventional result, 
Eq. (9). It is worth noticing that option B is Jn accordance with the quantum 
statistical mechanical result for the static mode a. 

In order to decide between the two options we have to bring physics into the 
consideration, i.e. the appropriate dispersion relation. This is the topic of the next 
subsection. 



3.2. The dispersive case 

Let d) be the frequency along the imaginary frequency axis. There are essentially 
two dispersion relations on the market. The first is the plasma relation 

S{iu;) = l + '^, (28) 

which is valid at the far ultraviolet for the light elements and at the X-ray region for 
heavier elements (Sec. 78 in Ref. 22). If we nevertheless employ Eq. (28) even near 
w = 0, it follows that n{iLo)uja/ 1^ — s- Xp/v where, in dimensional units, Xp = LOpa/c. 
Taking typically Wp ~ 3 x 10^^ s~^ and a ^ 1 cm we get Xp ~ 10^. In practice, 
the most significant values of I are much lower than this. We can thus assume that 
Xp/v S> 1, so that in practice 7 — > 0, 5 — > 0. That is, we recover in this way option 
A, and thereby the conventional result, Eq. (9). 

Consider next the Drude model for the dielectric, corresponding to 

e{^u) = l + —^^ (29) 

UJ{UJ + 7) 

7 being the relaxation frequency. According to this relation n{iil>)u) — s- when 
— !■ 0, implying that 7 — s- 1, 5 \. That is, we recover option B. The total to = 
free energy for a metal is thus according to the Drude model predicted to be one 
half of the conventional expression (9). 

Which of the two dispersion relations is correct? In our opinion, it is the rela- 
tion (29), when a) — > 0. On physical grounds the permittivity has to be inversely 
proportional to the frequency at low frequencies; cf. Sec. 77 in Ref. 22. Explicitly, 
£{iCj) = cr/w, where a is the conductivity. This is a result following directly from 
Maxwell's equations. The Drude model satifies this requirement. Thus both the 
Drude model (and, as we have seen, statistical mechanical methods), support the 
option B above. The plasma model, Eq.(28), as we have noted, is appropriate only 
at higher frequencies. 
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4. Conclusions, and Final Remarks 

We may summarize as follows: 

1. For a nondispersive dielectric, the formalism in spherical geometry becomes 
quite tractable. The free energy at finite temperatures is calculated from Eq. (2). 
Employing the Debye expansions up to 18th order in the quantity 6 (cf. Eq. (15)), 

good accuracy is achieved for all values of x and /. 

2. In the special case of a metal, adopting the Drude dispersion relation, we find 
that the m = TE mode does not contribute. The total m = free energy for a 
metal becomes accordingly only one one half of the conventional expression (9). 

3. As an extension of the above considerations, one may inquire about the 
magnitude of the m = contribution to the free energy for a dielectric. This case is 
treated in more detail in Ref. 3. From Eq. (2) one has, for arbitrary temperatures, 

/3F(m = 0) = - Y,{21 + 1) ln(l - X™) (30) 
;=i 

{Xf^ does not contribute for a dielectric). We may define Y as the ratio between 
F{m = 0) and the expression (2) for the full free energy: 

Y _ = 0) 
F 

For given d/a, Y thus becomes a function only of t. In Ref. 3, we show how Y 
varies with t for various values of d/a, for a fixed value of n. 

One general conclusion to be drawn from this calculation is that the less the 
value of d/a, the less becomes the importance of the m = term. This is a result 
that can be understood physically: when the slit is narrow, we can approximately 
regard the system as a conventional two-plates system. For the latter geometry, 
it is known that the classicality condition can be written as dT ^ 1, where d is 
the distance between the plates (cf. Sec. 82 in Ref. 23). When d decreases the 
system thus becomes more and more a quantum-mechanical system, necessitating 
an increasing large region of frequencies determining the value of F. The relative 
importance of the low frequencies, in particular that of m = 0, thus has to diminish, 
in accordance with the result of the calculation. 
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